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ABSTRACT 

We report deep Green Bank Telescope spectroscopy in the redshifted Hi 21cm and OH 18cm lines from the 
z = 0.765 absorption system towards PMN J0134— 0931. A comparison between the "satellite" OH 18cm 
line redshifts, or between the redshifts of the Hi 21cm and "main" OH 18cm lines, is sensitive to changes in 
different combinations of three fundamental constants, the fine structure constant a, the proton-electron mass 
ratio fi = m p /m e and the proton g-factor g p . We find that the satellite OH 18cm lines are not perfectly 
conjugate, with both different line shapes and stronger 1612 MHz absorption than 1720 MHz emission. This 
implies that the satellite lines of this absorber are not suitable to probe fundamental constant evolution. A 
comparison between the redshifts of the Hi 21cm and OH 18cm lines, via a multi-Gaussian fit, yields the 

strong constraint [AF/F] = [—5.2 ± 4.3] x 10~ 6 , where F = g p [/i« 2 ] 1 57 and the error budget includes 
contributions from both statistical and systematic errors. We thus find no evidence for a change in the constants 
between z = 0.765 and the present epoch. Incorporating the constraint [A/j//j] < 3.6 x 10~ 7 from another 
absorber at a similar redshift and assuming that fractional changes in g p are much smaller than those in a, we 
obtain [Act/ a] = ( — 1.7 ± 1.4) x 10~ 6 over a lookback time of 6.7 Gyrs. 
Subject headings: atomic processes — galaxies: high-redshift — quasars: absorption lines 



1. INTRODUCTION 

The standard model of particle physics implicitly assumes 
that the values of coupling constants and particle masses 
do not depend on space or time. Conversely, variation in 
such "fundamental constants" appears to be a generic feature 
of higher-dimensional theories ai ming to unify th e standard 
model and general relativity (e.g. Marciano 119841) . Studies 
of fundamental constant evolution are hence of much inter- 
est as they both probe the foundations of the standard model 
and allow the possibility of distinguishi ng between different 
unification models at low energy scales (Uzan 2 01 lb . 

While laboratory atomic clock studies have yielded strong 
constraints on short-term changes in the fine structure con- 
stant a (e.g. Rosenb and et al.ll2008l) . such studies are not sen- 
sitive to changes on Gyr timescales. A wide range of methods, 
based on various spectral transitions, has been used to probe 
changes in a, the proton-electron mass ratio /i = m p /m e 
and the proton g-factor on cosmological timescales 
(e.g.lBahcall et al.lll967t[T hompso n|[T97llWolfe et al.lll976t 
Varshalovich & Levshako vll 19931: iDzuba et al.lll999h |D arling 
2003 1 : IChe ngalur & Kanekar 2003; Kane kar & Che ngalur 
2004; Flambau m & Kozlovl 120071) . Significant progress has 
recently been made in both the development of new tech- 
niques and the sensitivity of measur ements. A co mbination of 
the many-multiplet method (Dzub a et al.ll 1999h with spectra 
from the High Resolution Echelle Spectrograph (HIRES) on 
the Keck telescope has found evidence for changes in a with 
redshift, [Aa/a] = (-5.7 ± 1.1) x 10~ 6 for 143 absorbers 
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at an average redshift z = 1.75 (Murphy et al. 2004). Later 
studies applying this method to spectra from the Ultraviolet 
Echelle Spectrograph (UVES) on the Very L arge Telescope 
(VLT ) have not confirm ed this res ult (e.g. iSrianand et al.l 
120071 iMolaroetal] 12008). Recently. IWebb et al.l ( fcOlHT ap^ 
plied the many-multiplet method to a large VLT-UVES sam- 
ple and also did not find evidence su pporting the Keck- 
HIRES result. While IWebb et al.l (120111) attempted to rec- 
oncile the Keck-HIRES and VLT-UVES results by propos- 
ing spatio-temporal changes in a, a simpler explanation is 
that the errors in both studies have been under-estimated, es- 
pecially given that systematic and unexplained errors have 
been shown to be p resent in the wavelength calibration of 
both spectrographs (iGriest et al.ll20ToHWhitmore et al.ll20iot 
lAgafonova et al.1201 11) . Keck-HIRES and VLT-UVES spectra 
in redshifted H2 lines have also yielded constraints on changes 
in n: the best cur rent result is [A/j/fi] < 4.4 x 10" 6 (z ~ 2.8; 
iKing et aT1l201 ll) . although the errors here too may have been 
under-estimated due to wavelength calibration issues. 

Radio spectroscopic techniques provide independent 
probes of fundamental constant evolution, with diff erent sys- 
temat ic effects from those in optical schemes (e.g. [Kanekar 
2008). Such methods include co mparisons between rota - 
tional and Hi 2 lcm hyperfine lines dDrinkwater et al.lll998l) . 
between differe nt hydroxy 1 (OH) lines or Hi 2 lcm and 
OH 18cm lines dDarljngl 120031: IChengalur & Kanekd |200"I 
iKanekar & Chengalurll2004l). between rotational and am mo- 
nia (NH3) inversion lines dFlambaum & Kozlovl |2007 | ). be- 
tween different methanol o r hydronium lin es (Jans erTet al.1 
1201 UlLevshakov et al.l20TltlKozlov & Levshakovl201 ll) . be- 
tween far-infrared fine structure lines and rotational lines 
dLevshakov et al.ll2008l) . etc, each sensitive to different com- 
binations of a, [i and g p . For example, the inversion-rotation 
comparison has yielded the best current constraint on changes 
in n from any astro nomical method , [A/i/n] < 3.5 x 10~ 7 
over < z < 0.685 (lKanekarll2011l) . 

In general, techniques using multiple spectral lines from 
a single atomic or molecular species (e.g., OH, methanol, 
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Fell, etc) are preferable to those using different species, as 
lines in the former case are likely to arise in the same gas, 
making the technique less susceptible to local velocity off- 
sets. The satellite OH 18cm lines have the special property of 
having exactly the same shape and opposite sign in certain 
astrophysical circumstances, due to a population in version 
mechanism and quantum mechan ical selection rules (Elitzur 
fl99l van Langeveld e et al.|[l995l) . This "conjugate" behav- 
ior ma kes them ideal probes of fundamen tal constant evo- 
lution dKanekar et alj|2004 iKanekari 120081) . Only two such 
conjugate OH 18cm systems have so far been discovered at 
cosmological distances, at z ~ 0.24 7 towards PKS 1413+135 
dKanekar et a lj |2004t lDarling|l2004l) and z ~ 0.765 towards 
PMN J0134-0931 dKanekar et al.ll2005l) . A high-sensitivity 
study of PKS 1413+135 with the Westerbork Synthesis Ra- 
dio Telescope and the Arecibo Telescope found tentative evi- 
den ce (at 99.1% confidenc e level) for changes in a, /i and/or 
g p dKanekar et al.l 1201 Oal) . In this Letter, we report deep 
Green Bank Telescope (GBT) observations of the redshifted 
Hi 21cm and OH 18cm lines in the z ~ 0.765 system towards 
PMN JO 134— 0931, that yield strong constraints on changes in 
the fundamental constants. 

2. OBSERVATIONS, DATA ANALYSIS AND SPECTRA 

The GBT observations of the Hi 2 lcm and OH 18cm 
lines from the z ~ 0.765 g ravita ti onal lens towards 
PMNJ0134-0931 dWinn et all 120021: iKanekar & Briggsl 
2003) were carried out between September 2005 and October 
2006 (proposal AGBT05C-037), using the PF1-800 and 
PF2 receivers, respectively. The observations used the 
Autocorrelation Spectrometer (ACS) with 9-level sampling, 
two polarizations and 10-second integrations. A single ACS 
12.5 MHz band, sub-divided into 32768 channels, was used 
for the Hi 21cm line, while four ACS 12.5 MHz sub-bands, 
each sub-divided into 8192 channels, were used for the 
ground-state OH 18cm lines. The system passband was 
calibrated by position-switching every five minutes, with 
system temperatures measured using a noise diode. The 
on-source times were ~ 15 hours for the Hi 2 lcm line and 
~ 30 hours for the OH 18cm lines. 

All data were analysed in the package DISH, using stan- 
dard procedures. The OH 18cm data were especially affected 
by intermittent radio-frequency interference (RFI); a visual 
inspection of every calibrated 10 second integration was used 
to excise data affected by RFI. Most of the data were found 
to have narrow-band RFI at ~ 945 MHz, close to the centre 
of the OH 1667 MHz absorption profile. The channels af- 
fected by RFI were entirely edited out in the final spectrum, 
along with 25 channels on either side (i.e. the frequency range 
944.867-945. 150 MHz). Similarly, weak RFI was also found 
adjacent to the satellite OH 18cm profiles, due to which a few 
channels in the final spectra were blanked out. 

The final Hi 21cm and "main" OH 18cm optical depth spec- 
tra are shown in the upper panels of Fig. [T] while the "satel- 
lite" OH 18cm spectra are shown in Fig. [2] The Hi 21cm and 
OH 18cm spectra have velocity resolutions of ~ 0.3 km/s and 
~ 1.0 km/s, respectively (after Hanning-smoothing and re- 
sampling), and root-mean-square optical depth noise values 
of 0.0019 per 0.3 km/s channel (Hi 21cm) and - 0.0011 - 
0.0012 per 1.0 km/s channel (OH 18cm lines). Note that 
these are "apparent" optical depths, derived using the to- 
tal flux density of PMN J0134-0931 at the line frequencies. 
PMN J0134— 0931 is a 5-component g ravitational lens sys- 
tem with an angular extent of re 0.7" dWinn et al.ll2.002h . (i.e. 



re 5kpc at z = 0.765 and unresolved by the GBT beam). 
The Hi 21cm and OH 18cm features are likely to arise against 
only one or two of the source components, implying that the 
"true" optical depths are probably significantly larger than the 
measured optical depths. 

Fig. |2] shows that, while the 1720 MHz line is in emis- 
sion and the 1612 MHz line is in absorption, the two lines 
have different strengths, with the peak optical depths in the 
1720 MHz and 1612 MHz lines being re -0.007 and re 0.01, 
respectively. Thus, although the satellite OH lines have simi- 
lar shapes, they are not exactly conjugate, with the 1612 MHz 
line about 1.5 times stronger than the 1720 MHz line. We tried 
a number of analysis procedures and RFI excision schemes to 
test whether the difference between the satellite line profiles 
might arise due to RFI; the difference was found to be present 
in all cases. Further, the flux densities of PMN J0134— 0931 
measured in the 1612 MHz and 1720 MHz spectra were very 
similar (~ 0.7 Jy in each), and the RMS noise values on the 
two spectra are comparable, indicating that there is no scaling 
error in the flux density calibration. It is thus unlikely that 
the difference between the satellite line profiles arises due to 
either RFI or problems with the flux density scale. 

3. PROBING FUNDAMENTAL CONSTANT EVOLUTION 
3.1. The satellite OH lines 

The 2n 3 / 2 (^ = 3/2) OH rotational ground state is split 
into four sub-levels by A-doubling and hyperfine splitting; 
two sub-levels have total angular momentum quantum num- 
ber F = 2, while the other two have F = 1. The satellite 
OH 18cm lines correspond to transitions with AF = ±1 be- 
tween these sub-levels. Similarly, two sub-levels of every ex- 
cited rotational state have F = J + 1/2, while the other two 
haveF = J - 1/2. 

The satellite OH 18cm lines are said to be "conjugate" when 
they have the same shap e, but with on e line in emission and 
the other in absorption (Elitzuri ri992l) . This arises because, 
when the OH molecules are pumped to excited states (by col- 
lisions or far-infrared radiation), the downward cascade to the 
ground state yields population inversion in the ground state 
sub-levels as certain transitions are forbidden by the selection 
rules AF = 0, ±1. If the last stage of the cascade is the intra- 
ladder 119/rni transition 2n 3/2 (J = 5/2) -> 2n 3/2 (J = 
3/2), transitions between the F = 3 and F = 1 sub-levels are 
forbidden; this would over-populate sub-levels with F = 2 
relative to those with F = 1. As a result, the 1720 MHz 
transition would be inverted while the 1612 MHz transition 
would be anti-inverted; this is the situation in Fig.lH with the 
1720 MHz and 1612 MHz transitions in emission and absorp- 
tion, respectively. 

If the 119/im lines connecting the excited and ground states 
are optically thick (which depend s on the local particle num- 
ber de nsity and velocity gradient; Elitzmi n"976t iGuibert et alJ 
1978), the rate coefficients for the different branches of the 
cascade are independent of line strength. This implies that 
the same number of particles are present in the two F = 2 
ground-state sub-levels, and, similarly, in the two F = 1 sub- 
levels. The 1720 MHz and 1612 MHz lines menhaye identi- 
cal strengths and shapes, albeit opposite sign (lElitzurl 11992). 
This is the ideal situation for the use of the satellite OH 18cm 
lines to probe changes in the fundamental constants as the 
identical l ine shapes guarantee that the lines arise fr om the 
same gas dKanekar & Che ngalur 20(3 |Kanikar 2008). 

Exactly conjugate satellite OH 18cm lines have been 
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FIG. 1. — GBT optical depth spectra in the redshifted [A] Hi 21cm and [B] "main" OH 18cm transitions; the solid line in each upper panel shows the 4- 
component fit, with the fit residuals shown in the lower panels. Shaded regions indicate frequencies affected by RFI. 

rise to part of the main OH 18cm absorption. For example, 
there is a main OH 18cm component at z = 0.76385 (see 
Table 1), with peak 1667 MHz opacity w 0.0129. If the 
OH ground-state levels in this cloud are thermalized, the 
satellite lines would have peak optical depths of m 0.0014 
(nine times weaker than the 1667 MHz line). The 1612 MHz 
absorption would hence be increased, and the 1720 MHz 
emission reduced, by this amount, yielding a difference of 
~ 0.0028 between the satellite optical depths, consistent with 
the observed difference. 

Since the satellite OH 18cm lines are not perfectly conju- 
gate, it cannot be assumed that the lines arise from the same 
gas. The lines also have slightly different shapes, in addition 
to the different strengths. This implies unknown systematic 
effects if these lines are used in studies of fundamental con- 
stant evolution. We hence conclude that the satellite OH 18cm 
lines in the z ~ 0.765 absorber towards PMN JO 134— 0931 
are not suitable to probe changes in a, /i and g p . 

3.2. The "main " OH 18cm and Hi 21cm lines 

A comparison between the redshifts of the "main" 
OH 18cm lines and the Hi 21cm line is sensitive to changes 

in F = g p [a 2 /i] ' (Cheng alur & Kanekarll2003l) . We used 
a multi-Gaussian fit to the Hi 21cm and main OH 18cm lines 
to test for changes in F, assuming the same velocity structure 
in the different lines. Independent fits to the Hi and OH lines 
found a 4-component model to give a good fit to each line, 
with good agreement between the line widths of the corre- 
sponding Hi and OH components. Turbulent broadening was 
hence assumed to dominate the line widths, with the Hi and 
OH widths of each component tied together in the fit. No as- 
sumption was made about the relative strengths of the main 
OH lines. However, since the main OH 18cm line frequen- 
cies have the same dependence on a, 11 an d g p to first order 
(IChengalur & Kanekar!l2003HKozlovll2009l) . the OH line red- 
shifts of each component were tied together. The fit included 
a single velocity offset between all Hi 21cm and OH 18cm 
components, to account for the putative change in F. 
With the above assumptions, we used VPFIT to carry out a 
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FIG. 2. — GBT optical depth spectra in the redshifted OH 1612 MHz line 
(upper panel) and OH 1720 MHz line (lower panel). Shaded regions indicate 
frequencies affected by RFI. 



observe d in several extragalact ic sourc es including 
CenA (Ivan Langevelde et all 11995b. M82 dSeaquist et all 
[1997b. NGC 253 dFraver et all 11998b. and PKS 1413+135 
dKanekar et al.ll2010al) . In Cen A and NGC 253, the satellite 
lines are conjugate over a wide range of conditions and even 
show the cross-over from absorption to emission (and vice- 
versa) in each transition. As such, the fact that the satellite 
lines towards PMN J0134— 0931 are not perfectly conjugate 
is both unusual and unexpected. If the 119/im transitions that 
dominate the downward cascade are not optically thick, the 
rate coefficients of the decay routes to the different sub-levels 
are different and one would not obtain conjugate behavior 
dElitzurlll976l) . Alternatively, there may be absorption at the 
satellite line velocities from the molecular cloud that gives 
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simultaneous multi-Gaussian fit to the Hi 21cm and OH 18cm 
profiles, varying the fit parameters to minimize \t- A four- 
component model yielded \t = 109 and noise-like residuals 
after subtracting out the fit. A Kolmogorov-Smirnov rank-1 
test found the residual spectra (after subtracting out the fit) to 
be consistent (within ~ 2a significance) with a normal dis- 
tribution. Fits with fewer than four components were found 
to yield a significantly higher xl', conversely, increasing the 
number of components did not improve the fit. 

A four-component model thus provides a good fit to the 
Hi 21cm and OH 18cm profiles. The number of fitted param- 
eters here is 21, four redshifts, four line widths, twelve peak 
line depths in the Hi 2 lcm, OH 1665 and OH 1667 spectra, 
and the velocity offset between the Hi 21cm and OH 18cm 
lines. The parameters of the fit are summarized in Table 1; 
the error on each parameter has been increased by a fac- 
tor yx?, to account for the fact that the best-fit value of 
xt is slightly larger than unity (x 2 , = 1.09), probably be- 
cause the RMS optical depth noise on the spectra has been 
marginally under-estimated. The best-fit velocity offset is 
AV^ = (-1.57 ± 0.44) km s _1 , with the Hi 21cm line blue- 
ward of the OH 18cm lines. 

The above error on AV is the statistical error from the 
fit. Other contributions to the error budget include errors 
in the line rest frequencies and the frequency scale calibra- 
tion, and local velocity offsets between the clouds giving rise 
to the Hi 2 lcm a nd OH 18cm lines; th e latter dominate the 
systematic errors. Kane kaFet alj (12005) estimated the Galac- 
tic dispersion between Hi 21cm and OH 18cm velocities to 
be < 1.2 km s -1 , using the measured dis persion between 
Galac tic Hi 21cm and HCO + velocities (Drinkwater et al. 
1998) and the good m atch between HCO + and OH velocities 
dLiszt & Lu cas 2000). We will assume that this dispersion 
also applies to the z ~ 0.765 absorber. For comparison, the 
main OH 18cm frequenci es have been measu red with an ac- 
curacy of 12 Hz (~ 2 m/s: lHudson et al.l2 006). while the error 
in the GBT frequency scale (mainly due to Doppler tracking) 
is < 15m/s. 

Our final result for the velocity offset between Hi 21cm and 
OH 18cm lines, including both statistical and systematic er- 
rors, is thus AV = [-1.57±0M(stat.)±1.2(syst.)]km&- 1 . 
This yields [AF/F] = [-5.2 ± l.5(stat.) ± A.0(syst.)] x 
10~ 6 . Adding the statistical and systematic errors in quadra- 
ture gives [AF/F] = [-5.2 ± 4.3] x 10~ 6 . We thus find no 
evidence for a change in a, /i or g p between z — 0.765 and 
today, i.e. over a period of 6.7 Gyrs0 

4. DISCUSSION 

The previous section constrains changes in F = 

g p [a 2 /i] 1 57 and does not yield independent constraints on 
changes in a, fj, and g p . However, this result can be used to 
obtain the sensitivity to changes in each parameter by assum- 
ing that the other two do not vary with time. The la sensi- 
tivities are [Aa/a] = 1.4 x 10~ 6 , [Afj,/fx] = 2.7 x 10~ 6 
and [Ag p /g p ] = 4.3 x 10~ 6 . A st ringent constra int on 
changes in [i was recently obtained by Kanekar (201 lj) from 
a comparison between NH3 inversion and CS/H2CO rota- 
tional lines at z = 0.685 towards B0218+357: [Afj,/fj] < 
3.6 x 10~ 7 between z = 0.685 and the present epoch. As- 
suming that there are no spatial variations in the constants, 

6 We use a standard LCDM c osmology, with Hp = 70.4 km s _1 Mpc -1 , 
n m = 0.27 and Q A = 0.73 (Komatsu et al. 201 1). 
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FIG. 3. — A comparison between the best present estimat es of [Aa/a], 
from a variety of methods: (I) conjugate satellite OH 18cm ( Kan ekar et al.1 
I2010d) . (II) Hi 21cm and OH 18cm lines (this work), (III) Hi 21cm 
and Cl lines jKanekar et al.1 I2010bt) . (IV) VLT-UV ES many-multiple t 
(Agafonova et al. 201 1), (V) VLT- UVES many-multipl et (Webb et alJ201 II) . 
(VI) Keck-HIRES many-multiplet (Murphy et al. 2004) and (VII) alkali dou- 
blet (Murphy etaLl BOOl). Methods (I), (II) and (III) assume [Afl/fj] « 
[Aa/a] and [Ag p /g p ] << [Aa/a]. Results from com parisons between 
the strongest or nearest met al and Hl-21cm components (e.g. Tzana varis et al.1 
I2005t ISrianand et al. 2010) hav e been exclud ed as they have unknown sys- 
tematic eiTors (e.g. Kanekar et al. 2006, 2007). 

we can replace [A/i//i] in the expression for [AF/F] to get 
[Ag P /g P ] + 3-14 x [Aa/a] = (-5.2 ± 4.3) x 10~ 6 . If 
we further assume that [Ag p /g p ] << [Aa/a], we obtain 
[Aa/a] = (-1.7 ± 1.4) x 10~ 6 between z = 0.765 and 
today. 

Fig. [3] shows a comparison between the bes t estimates of 
[Aa/ a ] from a variety of technique s t oday (|Murr>hv et all 
2001L I20041 iKanekar et al.1 l2010allbl: IWebbetalJ I20TTF 
Agafo nova et al.l 1201 ll) . The Keck-HIRES many-multiplet 
dataset is the only one that finds statistically-significant 
evidence for changes in a. However, systematic wave- 
length calibration e rrors may have been under-estimated here 
(Gries t"et al.ll2.010h . At present, there appears to be no strong 
evidence from astronomical spectroscopy for changes in the 
fundamental constants on cosmological timescales. 

In summary, we have carried out deep GBT spectroscopy 
in the redshifted Hi 21cm and OH 18cm transitions from the 
z = 0.765 absorber towards PMN J0134— 0931. We find 
that the satellite OH 18cm lines are not perfectly conjugate, 
with the 1612 MHz absorption ~ 1.5 times stronger than 
the 1720 MHz emission. The fact that the satellite lines 
have different shapes implies that these should not be used 
to probe fundamental constant evolution, due to the possi- 
bility of unknown systematic effects. We obtain tight con- 
straints on changes in the quantity F = g p [na 2 ] 157 via 
a simultaneous fit to the Hi 21cm and OH 18cm line pro- 
files. A four-component Gaussian model assuming turbulent 
line broadening yields a good fit to both profiles, yielding 
[AF/F] = [-5.2 ± 4.3] x 10~ 6 , including both statistical 
and systematic errors. We find no evidence for a change in a, 
[i or g p between z = 0.765 and today. 
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Velocity offset 
km/s 


Heliocentric redshift 


FWHM 
km/s 


( 

Hi 21cm 


Dptical depth X 10 
OH- 1667 


OH- 1665 


A 
B 
C 
D 


(-1.57 ±0.44) 


0.765057 (11) 
0.7648025(78) 
0.7638501(42) 
0.7635476(25) 


120.5 ±2.2 
47.7 ±4.3 
51.9 ± 1.6 

37.76 ±0.63 


1.795 ±0.055 
0.67 ±0.12 
1.064 ±0.043 
2.868 ±0.055 


0.509 ± 0.024 
0.321 ± 0.049 
1.291 ±0.038 
1.656 ±0.046 


0.325 ±0.019 
0.167 ±0.035 
0.754 ±0.028 
1.119 ±0.035 



TABLE 1 

The parameters of the best 4-Gaussianfit to the Hi 2 Icm and OH 18cm profiles; the second column contains the best-fit velocity 

OFFSET BETWEEN THE LINES. 
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